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ABSTRACT 

Narrow-line Seyfert 1 galaxies (NLSls) with very small broad-line widths (say, 
FWHM(H/3) < 1200kms -1 ) represent the extreme type of Seyfert 1 galaxies that have 
small black hole masses (M B h) and/or high Eddington ratios (L/L Edd ). Here we study 
the X-ray properties of a homogeneously and optically selected sample of 13 such 
objects, termed as very narrow line Seyfert 1 galaxies (VNLSls), using archival XMM- 
Newton data. It is found that the Fe Ka emission line is at most weak in these objects. 
A soft X-ray excess is ubiquitous, with the thermal temperatures falling within a strict 
range of 0. 1-0.2 keV. Our result highlights the puzzling independence of the thermal 
temperature by extending the relations to even smaller FWHM(H/3), i.e., smaller M BH 
(~ 10 6 M ) and/or higher L/L Edd . The excess emission can be modeled by a range 
of viable models, though the disk reflection and Comptonization models generally 
give somewhat better fits over the smeared absorption and the p-free models. At the 
Eddington ratios around unity and above, the X-ray spectral slopes in the 2-10 keV 
band are systematically flatter than the Risaliti et al.'s predictions of the relationship 
with L/L Edd suggested previously. Short timescale (1-2 hours) X-ray variability is 
common, which, together with the variability amplitude computed for some of the 
objects, are supportive of the scenario that NLSls are indeed AGN with relatively 
small M BH . 
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Introduction 



Type 1 active galactic nuclei (AGNs) are characterized by prominent broad emission lines in 
their optical/UV spectra. The lower end of the line widths is mostly populated by the so-called 
narrow-line Seyfert 1 galaxies (NLSls), defined as having the broad hydrogen Balmer lines nar- 
rower t±jm^2J)00kmjf^^ (FWHM) and the relatively weak [O III] 
lines (IQsterbrock & Poggel 19851 : lGoodrichlll989|) . NLSls show some extreme properties among 



Komossa 



2008 , for a recent review). Previous studies have revealed a set of corre- 



AGNs (see 

lations among AGN optical emission lin e and X-ray properties — the so-called eigenvector 1 (EV 
1) correlations (|Boroson & Greenlll992l) . which is believed to be driven by the Eddington ratio 
(L/L E dd). A narrow H/? l ine is generally associa ted with strong optical Fe II and weak \Q IIIl 



emiss i on (IGoo dric 



1996; IWang et al. 



n 



1989c IVeron-Cetty et al.ll2001l) . a steep soft X-ray spec tral slope (IBoller et al. 



1996 ). and fast and large amplitude X-ray variability ( Leighly 19991 : Grupe 



20041) . However, these correlations were found based on AGNs with FWHM(H/3) mostly broader 
than ~ 1000 km s" 1 , below which only a small number of objects were known by then. One 
would expect naively, based on the EV1 correlations, that NLSls with very small width (say, 
FWHM<1000kms _1 ) would show even extreme properties in X-ray, i.e., even steeper soft X-ray 
slopes and even faster and larger- amplitude variability. We refer to such AGNs as very narrow-line 
Seyfert 1 galaxies (VNLSls) hereafter in this paper. 

It has been recently found that the gas motion in the broad-line region (BLR) is virialized 



luminosity with an index of roughly 0.5 (IKaspi et al. 



([Peterson & Wandell l2000al : lOnken & Peterson! 12002) and that the BLR size scales with optical 



20051 : iBentz et al.ll2006l) . A combination of 



these properties naturally li nks FWHM(H/3 i ) with black hole mass M BH and L/L EM in a way as 
FWHM 4 ~ MBH^/^Edd)" 1 iMcHardvetalil2006l) . Therefore, a narrower FWHM(H/?) indicates 
a larger (L/L E dd)/^BH ratio, provided that the inclination is not a dominating eff ect. This is what 



NLSls are commonly thought to be, as argued ex tensively in the literatures (e.g jMineshige et al 



2000; IPeterson et al.ll2000bl : ISulentic et al.ll2000l) . and should be even more extreme for VNLSls 
as expected. 

As such, the X-ray properties of these extreme black hole accreting systems are of particular 
interest, in light of the following considerations. Firstly, VNLSls are well suited for investigating 
the soft X-ray excess emission commonly detected in Seyfert 1 galaxies and quasars, whose ori- 
gins remain controversial. One motivation is to attempt to link the soft X-ray excess with the black 
body emission from accretion disks, whose maximum temperature would be the highest among 
AGNs currently known (since r max oc[(L/L Edd )/M B H] 1/4 oc (FWHM)~ l ), and might be detectable 
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with the current X-ray satellites. Interestingly, this attempt was successful in at least one AGN, 
RX J1633 +4718, that is also a VNLS1 (FWHM(H/?)~ 900km s" 1 ) but radio-loud, as recently dis- 
covered by I Yuan et al.1 (I2010|) . though similar cases are extremely rare. Alternatively, the observed 
soft X-ray excess can be mimicked by relativistically blurred line emission of the refle ction compo 



nent f rom a highly ionized inner disk, which may be dominant in high L/L Edd systems (IFabian et al. 



20021) . Secondly, there were suggestions that NLSls resemble the fastest acc reting states ('high' 



and 'very high' states) o f X-ray binaries, in both the X-ray spectra (e. 



Middleton & Pond 120071) and X-ray quasi-periodic oscillations (QPOs) (|Gierliriski et al. 



Pounds et al. 



1995 



2008 



and thus VNLSls are more suitable for studying such an analogy. Thirdly, the X-ray spectral and 
temporal properties of V NLSls can be compa r ed to AGNs with genu i ne small mass black holes, 
say, M BH < 1O 6 M (e.g.. iGreene & Holl2007al : bewangan et al.ll2008l : iMiniutti et al1l2009h . This 



may provide possi ble insight in the black ho l e masses of VNLSls and help distinguish different 



models of NLS 1 s (lOsterbrock & Poggell 19851 : iMineshige et alJl2000b iSulentic et al.1 12000) 



h 



Although some VNLSls have been studied in X-rays individually in the literature, systematic 
studies of their ensemble X-ray properties are rare, however, given the lack of homogeneously 
selected samples in the past. Recently, the X-ray properties of small s amples of AGNs with 



Mjm< 10 6 M^with XMM-Newton observations have been presented (e.g.. iDewangan et al.1 12008 



Miniutti et al.l 120091) . among which several objects are in fact VNLS1 considering their optical 
spectral properties and the high L/L Ed d. It was found that these VNLSls are characterized by 
strong and rapid X-ray variability and soft X-ray excess emission. However, more observations for 
a larger, homogeneously selected sample are needed to confirm these results. 



Using a large NLS1 sample selected from the SDSS. IZhou et al.l (120061) found that, to one's 
surprise, the previously known FWHM(H/3)-T s (soft X-ray photon index) anti-correlation becomes 
flattened at FWHM ~ 1200kms _1 . Though in the small FWHM regime AGNs having flat T s 
have previously been noted to exis{l the lack of expected steep soft X-ray slopes is intriguing. 
However, enhanced X-ray absorption in VNLSls may expla in such a trend, which, though seems 
to be unlikely, cannot be ruled, since in IZhou et al.1 (|2006l) the T s (estimated from the ROSAT 
hardness ratios) are subject to large uncertainties. Detailed spectral modeling of the X-ray spectra 
with higher spectral resolution and S/N for these objects is needed to confirm this interesting trend. 

Motivated by the above considerations, here we present a study on the X-ray properties of a 
sample of NLS1 with extreme linewidth, FWHM < 1 200km s" 1 , using data from archival XMM- 
Newton observations. The sample and data reduction are described in Section 2. The modeling of 
the XMM-Newton spectra are presented in Section 3, with a focus on the soft X-ray excess. The 



'For such A GNs, their deviation from the above relation is explained as due to their low luminosity (low L/LEdd); 
see iLaoj hood) . 
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X-ray temperal properties are presented in Section 4, followed by the effective optical to X-ray 
spectral indices. The results and their implications are discussed in Section 6, and are summarized 
in Section 7. We assume a cosmology with H = 70km s" 1 Mpc" 1 , = 0.73, and Vl M = 0.27 '. All 
quoted errors correspond to the 90% confidence level unless specified otherwise. 



2. Sample and X-ray data 



2.1. X-ray VNLS1 sample 



We select VNLSls from a large, homogeneous sample of ~2000 NLSls built by IZhou et al. 



(120061 hereafter Zhou'06 sample) from the SDSS DR3, which can be considered as basically 
optically selected. We adopt an operational linewidth cutoff of FWHM(H/3) < 1200 km s" 1 for 
VNLSls, in consideration of the fa ct that below roug hly this value the r s -FWHM anti-correlation 
becomes flattened (see figure 17 in lZhou et al.ll2006l) . There are 384 NLSls mee ting this criterion. 
We match these VNLSls with the 2XMM source catalogue (IWatson et al.ll2007l) using a matching 
radius of 5" and then select those detected in X-rays with at least 200 net source counts. We 
consider radio-quiej^l objects o nly, since X-rays from radio-loud NLSls may be con taminated by 
emission from relativistic jets (IZhou et al .1120071 : 1 Yuan et al.l 12008c lAbdo et al J 120091) . The above 
selection results in 13 objects with reasonable signal-to-noise (S/N) ratios, which form our working 
sample of VNLSls in this study. The sample objects are listed in Tabled] and the logs of the 
XMM-Newton observations are summarized in Tabled Among the sample, the XMM-Newton 
data of seven objects are presented here for the first time; while the XMM-Newton spectra of six 
objects^ have been presented previously in various details in the literatures for different aims. For 
the purpose of sample study using homogeneously derived results, we also re-analysis the XMM- 
Newton spectra of these objects, in the same way as for the other objects whose XMM-Newton data 
are presented for the first time here. 



T he optical spectral and continuum parameters of the sample objects are taken from lZhou et al. 
(|2006|) and given in Tabled] • The black hole masses are esti mated from the broad c omponent 



of the Ha line using the M B H-hne width-luminosity relation in iGreene & Hoi (I2007bl). We also 



estim ate the Eddington ratio L/L E dd assuming the bolometric luminosity as 9AL 5 i o (|Elvis et al. 



1994b . where AL 510 o is the monochromatic luminosity at 5 100 A. Figured] shows the distributions 



2 Havin g radio-loudness l ess than 10, defined as the rest frame flux ratio between the radio 1 .4 GHz and the optical 
g-band (see lZhou et al . 2006) 

3 They are J010 7+1408, Jl 14 0+0307, J1357+65 25 jPewangan etail l2008t iMiniutti et al.1 Eooj) . J1246+022 
dPorauet et al.ll2004 . J2219+1207 dGallo et alJ2006bh , and J1415-0030 dFoschini et alJ2004 . 
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of several parameters of our working sample, namely, M B h> L/L E dd and the optical Fe II emission 
multiplets strength i?457(E, in comparison with those of the NLS1 sample of Zhou'06 as well as 
the overall FWHM < 1200 km s" 1 (VNLS1) sub-sample. A few remarks can be made concerning 
the bulk properties of the sample. Firstly, both of these two VNLS1 sub-samples have similar 
i?4570 distributions to that of the parent sample, confirming their typical NLS1 nature. The general 
strong Fe II emission is also demonstrated in the composite SDSS spectrum of our XMM-Newton 
VNLS 1 sample (Figure©. Secondly, our VNLS 1 s have lower M BH and higher L/ L Edd distributions 
in general than the overall NLS1 sample, as expected from their narrower linewidths. Thirdly, our 
XMM-Newton sample is roughly consistent with the overall VNLS Is in these distributions. Thus 
our X-ray VNLS1 sample is not biased from, but rather representative of, optically selected NLSls 
with the smallest linewidth. This should be kept in mind when comparing our results with those 
obtained in previous studies, especially for X-ray selected NLS1 samples. 



2.2. X-ray observations and data reduction 

The observational data with XMM-Newton were retrieved from the XMM-Newton science 
archival center. For all but one object the observations were operated in the full window mode. For 
the only exception, J2219+1207, the MOS cameras were operated in the small window mode, in 
which a considerable fraction of the source counts in the wing of the point spread function (PSF) 
were lost, and thus only the PN data are used. The PN observation of J 1246+0222 experienced a 
pile-up, which is corrected by excising the core of the PSF with a radius of 10". Some of the data 
of the individual cameras as listed in Table[2] cannot be utilized due to the sources being either at 
the edges or in the gaps of the CCDs, out of the field of view, or on a bad CCD column. 

For XMM-Newton data reduction we use the standard Science Analysis System (SAS, v8.1.0.). 
The Observation Data Files (ODF) are processed to create calibrated events files with 'bad' (e.g., 
'hot', 'dead', 'flickering') pixels removed. The time intervals of high flaring backgrounds con- 
tamination are identified and subsequently removed following the standard SAS procedures and 
thresholds. Source counts are extracted from a circle with a radius ranging from, depending on the 
source position on the detector, 30" to 65" at the source position, and the background counts from 
a source-free region with a usually larger radius. To extract X-ray spectra only X-ray events with 
the pattern < 4 for PN and < 12 for MOS are used. Background subtracted light curves are also 
extracted from cleaned events files and are subsequently corrected for instrumental effects (such 



4 Defined as the Fe II(AA4434-4684) to H/3 flux ratio, where Fe II(AA4434-4684) denotes the flux of the Fe II 
multiplets integra ted over the wavelength range of 4434^4684 A, and H/3 the flux of the broad component of H/3; see 
Zhou et all (l2006l) . 
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as vignetting and dead time) using the S AS task 'epiclccorr' . 

When both available the two MOS spectra of each observation are all found to be consistent 
well with each other, and are therefore co-added to form a single MOS spectrum to increase the 
S/N (using the FTOOLS addspec 1.3.0). We set the low-energy cutoff of the s pectra to 0.2 keV , 
following the recommendation of the most recent EPIC calibration status report (IGuainazzill2009|) . 
The high-energy cutoff is set as such above roughly this energy the background spectrum starts to 
dominate. Since there is a diverse range of the spectral S/N, the high-energy cutoff varies among 
the sample objects. For four objects with low spectral S/N, only the 0.2-2.4 keV range is used. For 
the remaining nine objects, the high-energy cutoff of either 7 keV or lOkeV is used, depending on 
the spectral quality (see Table©. The EPIC spectra are then grouped in a way that t here are at leas t 



25 counts in each energy bin. Spectral fitting is performed using XSPEC (v.l2.3. lArnaudl ll996). 
Whenever both PN and MOS spectra are available for an observation, we perform joint spectral 
fitting with all the same spectral parameters tied together except the normalization. 



3. X-ray spectra analysis 

3.1. Continuum shape and Fe Kq emission line 

In order to compare the X-ray continuum slopes with results from other previous AGN studies, 
we first characterize the X-ray spectral shape with an absorbed power-law model in both the soft 
(0.2-2.4 keV) and "hard" (2-10 keV) bandfl respectively. The results are listed in Table[3j In the 
soft X-ray band, the model with a neutral absorption column density (iV H ) fixed at the Galactic 
value (A^ al ) yields acceptable fit for about half of the sample objects. In the remaining objects 
the fit can be improved by adding an extra neutral absorber in the objects' rest frame. The fitted 
excess absorption 7V H are small, however, comparable to N^ al . In only one object, J1415-0030, 
ionized absorption is required to yield acceptable fit, with an edge-like feature around 0.6 keV (see 
Section |3T2l) . We conclude that intrinsic absorption is not significant in these objects. We thus 
suggest that the observed flatte ning of the FWHM(H/3)-T s relation below FWHM ~ 1200 km s" 1 



as found in iZhou et al.l (120061) . is not caused by X-ray absorption, but most likely an intrinsic 
property. For those having more than one measurement, there seems to be little or no changes in the 
soft X-ray spectral shape, and the mean T s are calculated. The fitted T s values range from 2.03 to 
3.72. We quantify the intrinsic distribution of T s (assumed to be Gaussian) that is d isentangled from 



measurement errors using the maximum-likelihood method as first applied by iMaccacaro et al 



5 For three objects only the 2-7 keV band is used since the spectra above 7 keV are dominated by backgrounds; see 
Tabled 
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(119881) : we find a mean (r s ) = 2.83^, and a standard deviation a = 0.3 l+g fg (90% confidence), 



whose confidence contours are shown in Figured 

In the hard X-ray band there are nine objects having high enough spectral S/N ratios for 
measuring photon indices T h . The absorption Af H is fixed at the Galactic value. A power-law 
can well reproduce the observed hard X-ray spectrum in five objects, whereas the remaining four 
objects (J0922+5120, Jl 140+0307, J1246+0222, and J2219+1207) show a possible broad excess 
emission feature in the residuals around 5 keV. The fitted T h values are in the range of 1.95-2.39. 
The maximum-likelihood intrinsic distribution of T h has a mean of (r h ) = 2. 19^ }g, which is flat for 
typical NLS Is, and a small intrinsic scatter, a = O.O+q jf (see Figure[3]for their confidence contours). 

No Fe Ka emission line feature appears to be present in all the objects except J1357+6525, 
which shows a marginal feature of a narrow-line at around 6.4 keV. For this object adding a Gaus- 
sian (to an absorbed power-law model in the 2-10 keV band) improved slightly the fit, though with 
only A% 2 =5 for 3 degrees of freedom (dof), i.e. the addition of an extra Gaussian component is 
not statistically significant. Thus only upper limits can be derived on the equivalent width (EW) 
of any potentially narrow Fe Ka line at 6.4 keV (assuming a = 10 eV). The derived upper limits at 
the 90% significance are given in Table[3j We c ompared our obje cts with the Fe Ka line EW and 



X-ray luminosity relation for AGNs as given in lPage et al.1 (120041) . and found that the derived line 
EW limits are well consistent with the relation. 



3.2. Soft X-ray excess 

A comparison of the soft and hard X-ray spectral indices obtained above indicates an overall 
spectral steepening toward low energies in all of the objects, suggesting the presence of the soft 
X-ray excess. As demonstration, we show in Figure|4]the X-ray spectra of the 4 VNLSls which are 
presented for the first time, and the extrapolation down to 0.2 keV of the power-law model fitted 
in the hard X-ray band. Significant excess emission in the soft X-ray ban d is prominent, similar to 



that reported in the other objects of the sample (e.g.. lMiniutti et al.ll2009h . which is also confirmed 
here. We thus conclude that the apparent soft X-ray excess emission is ubiquitous in our VNLS1 
sample. 

Given the small M B h and high accretion rate (L/L E dd) in VNLSls, the expected blackbody 
emission from accretion disks is shifted toward higher energies compared to classical AGNs with 
more massive black holes, and the high energy turnover may start to emerge in the soft X-rays 
(e.g., kT„ Mr ~72eV for J0940+0324 assuming a Schwarzschild black hole, e.g. Peterson 1997). 
Thus the blackbody emission directly from the disks might be detected. We first model the soft 
X-ray excess with a blackbody. The model yields acceptable or marginally acceptable fits for all 
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except for two objects. For J0107+1408 an additional neutral absorber is required to improve the 
fit. For J 141 5-0030 a moderately ionized absorber is needed; adding an absorption edge improves 
the fit significantly with Ax 2 =12 for 3 additional free parameters. The fitted edge energy is 
0.67 ±0.03 keV with an optical depth of 0.55±0.2, corresponding to the K-shell binding energy of 
ionized Oxygen ions. 

The fitted temperatures are in the range of 1 00-200 e V (TablelH), in agreement with those 
found in AGNs having more massive black holes of 10 7 " 9 M jGierlMski & Donell2004l: iPorquet et al.1 



20041 : ICrummy et al.l 120061: iBianchi et alj|2009h . Hence our results confirm the extensi on of the 
canon ical 1 00-200 eV temperatures down to AGNs with M B h as low as around ~ 1 6 M Q (|Miniutti et al 
2009|) by adding more objects in this M B h range. The result is clearly demonstrated in Figure[5l in 
which our results are c ompared with those of AGNs and quasars with more massive black holes 
(|Piconcelli et al.ll2005l : ICrummy et al.ll2006|) . These values are still systematically higher than the 
maximum temperatures predicted for standard accretion disks. The independence of the thermal 
temperature on M B u over such a wide M B h range argues against the direct blackbody emission 
from accretion disks as the origin of the obser ved soft X-ray excess for the vast majority of AGNs, 
except for RX Jl 633+47 18 dYuan et allboiol) . 



Recent studies suggest that, similar to the blackbody temperature, the relative strength of the 



soft X-ray excess al so falls within a relatively small range for PG quasars (IPiconcelli et al 



and Seyfert 1 AGNs ( Middleton & Donell2007 ) and AGNs with small masses ( Miniutti et al 



2005) 



2009) 



Here we examine this quantity for the VNLSls of our sample. We estimate the relative strength as 
the luminosity ratio of the excess component, modeled as a blackbody, to the total luminosity in the 
0.5-2 keV band. The values are in the range of 8 %-3 8% with a mean of 21% (TablelU). Apparently, 
when combined with previous results, as shown in Figure[6l there seems no strong dependence on 
H/3 linewidth over a large range, FWHM(H/3) = 600-10,000 km s" 1 (a Spearman correlation test 
probability of 0.79). The mean relative strength of our VNLSls is somewhat smaller than that 
(31%) of the PG quasars in (|Piconcelli et al.ll2005l) derived from their fitting results, though further 
confirmation is needed given the relatively small size of our sample. 

There are currently several viable models to account fo r the soft X-ray excess. P hoton trap- 
ping in high accreting system where advection is important (|Abramowicz et al.ll 19881) . or Comp- 
tonization of ultraviolet photon from the accretion disc by electrons as hotte r skin above the dis c 
dCzernv & Elvisll 19871: Iwandel & Petrosianll 19881 Ishimura & Takaharalfl993l : Iczernv et alJl2003h . 
can explain the required higher temperature. On the other hand, absorption or emission lines 
arising from atomic processes, when blurred due to relativistic motion, can mimic the soft X-ray 
excess. For example, strong relativistically-blurred emission/absorption lines between ~0.7-2 keV 
due to OVIFOVIII and Fe transition can be produced from ionized disc illuminat ed by an under- 
lying hard X-ray continuum (reflection) in the vicinity of the central black hole (|Ross & Fabian 



-9- 



1993|) . or from disc winds (|Gierlinski & Dondl2004l) . Below we investigate these models by fitting 



them to the soft X-ray excess spectra. In the spectral fitting, a power-law continuum modified 
by neutral absorption with A^ H fixed at the Galactic value is always included. For objects whose 
data have been analyzed previously, we compare our results with previous results individually in 
Appendix A. 



3.2.1. Comptonization 



We use the Comptonization model (comptt in XSPEC, lTitarchuklll994|) and fix the input seed 
photon energy at the innermost temperature of the standard accretion disk based on the estimation 
of the black hole masses and accretion rates (L/L Edd ). In this case the emergent spectral shape 
depends on only two parameters, the temperature and the optical depth of scattering electrons. 
This model gives significantly improved fits over, or at least as good as, the above blackbody fits 
for all of the objects (see TableH Figure[7]). The inferred electron temperatures are found in a 
relatively small range, fcr plasma ~ 0.17-0.30keV, and the optical depth, r ~ 10-25. 



3.2.2. disk reflection 



We use the latest ionized disk reflection model from lRoss & Fabianl (|2005|) {reflion in XSPEC) 
and in the fits the photon indices of the ionizing continuum and the observed continuum are tied to- 
gether, and the solar ab undance is assumed. For relativistic blurring the laor kernel model (kdblur 



in XS PEC . lLaorl 119911) is used with an outer radius fixed at 400 r g , an emissivity index of the disk 
fixed to the standard value of 3, and the inner disk radius allowed to vary. Although this model 
can reproduce acceptable fitting results for most of the objects, in three objects t he residuals in 



the so ft X-rays indicate possible contribution from another component. Following iMiniutti et al 



(|2009|) . we then include an additional blackbody component in the fits to account for possible con- 
tribution from the accretion disk. This improves the fits for all of the three, namely JO 107+ 1408 
(A X 2 /dof = 14/3), J0922+5120 (A X 2 /dof = 204/3), and J2219+1207 (A X 2 /dof = 24/2), with 
the addition of the blackbody component is statistically significant (a probability level <0.05 us- 
ing the F-test). The inferred temperatures are 0.05±0.008keV, 0.06±0.02keV, and 0.08^ keV, 
respectively, broadly consistent with the maximum temperatures at near the inner disk predicted 
from the estimated black hole masses and accretion rates. 

The disk reflection model, either with or without additional blackbody emission, provides 
acceptable fits for all and the best fits for some of the objects (see Figure[7]and Table|4]). The best- 
fit disk inner radius is less than 6r g , suggesting a highly spinning Kerr black hole for most of our 



objects. The inferred disk inclination varies from 0° to 50°. The ionization parameters are log£ ~ 
3.13-3.83. We define a parameter 'flux fraction' as the fl ux ratio of the reflec ted component to 
the total component over the range 0.2-10 keV (following ICrummy et al.ll2006l) . As can be seen 
in TableHl in most of the objects where this model gives a good fit, the reflection component 
dominates largely the total observed flux in the 0.2-10 keV band. 



3.2.3. smeared absorption model 

Finally we fit the spectra with the relativistically smeared absorption model (swindl in XSPEC). 
This model provides acceptable fits for nearly half of the objects but not for the remaining objects 
(TablelU Figure©. The inferred column densities are in the range of _/V H ~ 0.9-5 x 10 23 cm" 2 , the 
ionization parameters log£ ~ 2.75-3.59, and the smearing terminal velocities are very high, close 
to 0.5 c. 



3.2.4. p-free model 



In this work we also try the p-free disk model (diskpbb in XSPEC) to account for the soft 
excess. For some of the objects, this model gives as good fits as the Comptonization model, with 
the inferred temperatures at the inner disk radius of 0. 1 5-0.35 keV, and the index of the temperature 
profile p = 0.32-0.64. However, the innermost disk radii derived from the fitted normalization are 
significantly less than the estimated gravitat ional radii, even af ter the correction for the spectral 
hardening factor (~ 1.9 as recommended by Kawaguchil (|2003l) for AGNs with Eddington ratios 
around 1). We thus consider the p-free model to be physically unrealistic and hence no fitted 
parameters are listed here. 



3.2.5. Summary of soft X-ray excess modeling 

In general, the soft X-ray excess of these VNLS 1 s can be reproduced by more than one model, 
which often cannot be distinguished based on fitting statistics. However, disk reflection and/or 
Comptonization are much more preferred than the other two ones, i.e., smeared absorption model 
with marginally improved fitting and p-free model which is physically unrealistic. In the mod- 
eling the X-ray spectra, no additional intrinsic absorption are required for all the objects except 
J0107+1408 and J1415-0030, of which additional neural and warm absorption needed, respec- 
tively. 
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For SDSS J0107+1408, neutral absorption with N H ~6.55-1 1.41 x 10 20 cm" 2 is required. For 
J141 5-0030, in the above spectral fittings for this object with various soft X-ray excess models, we 
add ionized absorption {zxipcf'm XSPEC) applicable to all the emission components. The overall 
spectrum and the edge feature can be well reproduced and all the fits are improved significantly, 
with a decrease of A% 2 = 10 in general. The covering fraction of the absorber is close to unity, the 
absorbing column density is in the range of 9.3-15.8 x 10 20 cm" 2 and the ionization parameter of 
100.21-0.38^ depending on the exact model for the soft X-ray excess. 



3.3. IVFWH1VI relation 



We have showed above that there is little or no significant intrinsic abso rption in the X-ray 
spectra of most of these VNLSls. Therefore the photon indices derived in IZhou et al.l (I2006Q 
from the ROSAT ha rdness ratios are mostly reliable, and hence the spectral flattening towards the 
lower F WHM end (IZhou et al.ll2006l) is likely real, rather than being caused by X-ray absorption. 
Recently. iGrupe et al.1 (120101) studied the spectral indices of a sample of soft X-ray selected AGNs 
measured with the Swift XRT, some of which also have linewidt hs similar to ours. W e compare 
the soft X-ray photon indices of our VNLSls with those of the IGrupe et al.1 (|2010h sample, as 
shown in Figured It can be s een that in th e lowe st linewidth regime our T s values are statistically 
compatible with the result of Grupe et al. ( 201p[) @. It also appears that there is a lack of AGNs 
having both narrow FWHM (< 1000 km s" 1 ) and very steep soft X-ray spectra (T s > 3.5). This 
is consistent with the result suggested by IZhou et al. (bood) . Using the combined data points in 
Figured! we test explicitly whether there exists a significant flattening in the r s -FWHM relation 
at FWHM< 1000 km s" 1 using various methods; the result is inconclusive in the statistical sense, 
however. This might be partly due to the relative small sample size and/or the heterogeneity of 
the combined samples. A larger and homogeneously selected sample is needed to test the possible 
deviation of this well known r s -FWHM relation at the low-FWHM end in the future. 



Note that the T s values in IGrupe et al. are measured in the 0.2-2 keV band, slightly different from 0.2- 

2.4 keV adopted here; however, we find from spectral fits that the differences in r s thus caused are negligible for our 
objects. 
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4. X-ray variability 



Figure[9] shows, as examples, the 0.2-10keV lightcurves for the five objects in our sample, 
which are presented for the first time0. In fact, we find that, for most of the sample objects, the 
X-ray flux varied by more than a factor of 2 on timescales of 1-2 hours within the observational 
intervals. We conclude that short timescale variability is common to VNLSls. Remarkable flux 
variations in short-timescales are evident. We also investigate possible spectral variability within 
the observational intervals, which are divided into time bins, using the hardness ratios; however, no 
conclusive remarks can be made mainly due to the insufficient S/N of the data for such a purpose. 

The X-ray variability amplitude can be characterized by the excess variance^, which has 
been found to be str ongly correlated with the black hole masses (|Lu & Yull200ll : |Papadakisll2004l : 



O'Neill et al.l 120051). Using the XMM-Newton data, some are included in the current sample, 
Miniutti et all (|2009l) extended this relation to AGNs with black hole masses < lO 6 A/ (fhree in- 



cluded in our sample) and demonstrated that the relationship is relatively tight. This result indicates 
that black hole mass is a primary para meter that drives the relative X-ray variability in AGNs. In the 
same way as in lMiniutti et al.1 (120091) . we calculate the excess variance for the four objects not pre- 
sented previously, that makes use of lightcurve segments of equal duration (20 ks) and equal time 
bin size (500 s) in the 2-10keV band. However, for only one object, J0922+5120, the lightcurve 
has S/N high enough to yield reliably determined excess variance, logcr^xs = —1.15 ±0.68. We 
locate this obj ect (M RH = 10 6 63 M C .{ ) on the excess variance vs. black hole mass relation presented 
in Figure 8 of iMiniutti et al.1 (120090 . and find that it does follow closely th e relation. It should be 
noted that, Jl 140+0307, one of the three objects in IMiniutti et al.1 (12009b that are in our sample 
is also typical NLS1. As is generally believed that M BH is the underlying physical parameter that 
drives the dependence of the X-ray variability, the fact that these NLSls follow the same excess 
variance-M BH relation as for normal broad-line Seyfert 1 galaxies (BLSls) and quasars tends to 
validate their M BH estimation. That is to say, the black hole masses of the VNLSls in our sample 
are indeed relatively small, i.e. their narrow Balmer linewidth is caused primarily by relatively 
small M BH , rather than by a face-on flattened BLR as claimed in some papers in the literatures. 



^ForJ0107+1408 Jl 140 -0307, andJl 357+6525 t he X-ray lightcurves h ave been presented in lMiniutti etaD J2009h 
and bewangan et alj j2QQ8h : for J1415-0030 it was in lFoschini et alj J2004 . 

8 Defined as variance of the lightcurve in a specific time series after subtracting the contribution expected from 
measurement errors (e.g., Vaughan et al.lll999 ). 



-13- 



5. Optical/UV to X-ray spectral index 



The broad-band spectral energy distribution of A GN is commonly param eterized by the a ox 
parameter, a ox =-0.3838 log[L I/ (2500A)/L I/ (2keV)] (|Tananbaum et al.lll979). which is claimed 



to be correlated with the optical-UV luminosity fo r Seyferts and quasars (|Vignali et al.1 12003 



Yuan et al.ll 19981 : IStrateva et al.ll2005h . iGallol (12006al) reported that the objects in their NLS 1 sam- 
ple, which have systematically broader linewidths than ours, also follow the same relation for 
BLSls. Here we investigate this relation for our VNLSls. The flux density at 2500 A is calculated 
from the SDSS u band (effective observed-frame wavelength of 3543 A) PSF-magnitude adopt- 
ing the spectral slope of the composite SDSS quasar spectrum (ct„ = -0.44; Vanden Berk et al. 
2001). The a ox values of our objects are found to lie in the range from -1.41 to -1.18. We show 
in FiguredO] the a ox vs. 25 00 A monochr omatic luminosity relation for our objects, along with 
the normal NLS Is fro m the IGallol (I2006ah sample, as well as the regression relation for normal 



Seyferts and quasars (|Strateva et al.1 120051) . It shows that the VNLSls do follow closely the a ox - 
L uv relation defined by BLSls, and consistent with the result for NLS Is with broader linewidth 
(FWHM > 1000 km s" 1 ). 



6. Discussion 



6.1. Soft X-ray excess in VNLSls 



The ubiquity of the soft X-ray excess in our VNLS1 sample is interesting, which may have 
something to do with the generally high L/L E dd in our sample. The reflection model has been 
found to be a goo d description of the complexity of the X-ray spectra and spectral variability 



for se veral NLS Is (IMiniutti & Fabianll2004l : ICrummv et al.ll2006l : Ballo et al.ll2008l : IZoghbi et al 



2008J). In some cases the reflection component is found to dominate the observed X-ray band, 
which can be explained ei t her in terms of a corr ugated disk or strong gravitational light bending 
effects (|Fabian et al]|2002l) . IMiniutti et all (|2009l) found that this model reproduce well the XMM- 
Newton spectra of several AGNs with M B h < 1 6 M q , three of which are included in this paper. The 
spectral fitting for the other objects in our sample also supports this result. The small inner radii 
of the accretion disk derived argue for fast rotating black holes in most of these VNLSls, which 
may be a consequence of their fast accretion process. In a few cases the disk thermal emission is 
required, which is not unexpected given the relatively high disk temperatures for small M B h; this 
indicates that the disk reflection model is self-consistent. 



The Comptonization model successful explains t he soft excess in norma l Seyferts (IGierlihski & Done 



20041) and the spectral variability in RE J 1034+396 (IMiddleton et al J 12009) . However, as pointed 
out by IGierlihski & Pond (|2004l) . the derived electron temperatures and the optical depth are both 
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found in a small range (kT e ~ 0.1 -0.3keV and r ~ 10-20), which is puzzling and requires fine 
tuning of the disk parameters. Our result confirms that this is also the case for VNLSls. 

Although the 'p-free' model gives statistically acceptable fitting results in most of the objects, 
the inferred radii of the innermost disks are unphysically smaller than the gravitational radii. We 
thus consider the simple 'p-free' model to be unfeasible. The smeared absorption model can also 
reproduce the observed soft X-ray excess to some extent; however, it is disfavored since it results 
in generally poorer fitting statistics compared to the other models for our sample. Furthermore, the 
derived velocities are largely unconstrained at the extreme value and the ionization parameters are 
always found to lie in a narrow range. 



6.2. Correlations of spectral indices 

Our result is supportive of the lack of objects with very steep soft X-ray slopes and FWHM(H/3)< 
1200 km s" 1 in the r s -FWHM diagram; those objects may be expected simply based on this well- 
known EV1 correlation. Here we try to link this relationship with L/L Edd which is believed to 
be the underlying driver of the EV1 correlations. The soft X-ray spectral slope is primarily deter- 
mined by two factors, the slope of the underlying (hard X-ray) continuum and the effect of the soft 
X-ray excess. We examine the effect of the former by comparing the soft and hard band spectral 
indices of our sample objects, and find that they are strongly correlated with each other (a Spear- 
man correlation test probability P < 10" 5 ). Thus we conclude that soft X-ray slope T s is largely 
determined by, or tracing the underlying continuum slope T h . 

Compared to the T-FWHM relation, a more significant, and perhaps fundame ntal, correlation 



is the Th-L/LpM correlati on that has been detected in the range of L/L Edd < 1 i.e., (IShemmer et al. 



20061 : iRisaliti et alJl2009h . We show in Figured!] the r h -L/L Edd rela tion for our VNLSl s whose 



Th are available. Also plotted is, for a comparison, the regression of IRisaliti et al.l (120091) for the 
strongest Fh-L/L Edd relation (M B h estimated from H/3, the same as in our paper) and its extrapola- 
tion to the high L/L Edd regime where our sample objects are located. Interestingly, the T h values of 
all our VNLSls fall systematically below the extrapolation of the Risaliti et al.'s relation to high- 
L/L Edd values. To check whether this flattening might be caused by the enhanced Compton reflec- 
tion hump — known to exist in Seyfert galaxies and to make the hard X-ray spectrum flatten — in 
high-L/L Edd objects, we over-plot the 'underlying' X-ray continuum slope inferred from the above 
disk reflection model fitting. As can be seen, the 'underlying' hard X-ray continua (open circles) 
are indeed slightly steeper than the 'observed' one, but are still systematically flatter than the pre- 
diction of the Risaliti et al.'s relation. We thus suggest that, at L/L Edd ~ 1 or above, the hard X-ray 
spectra become flattened than what the Risaliti et al.'s relation predicts, at least for the VNLSls in 
our sample. Certainly more observations are needed to confirm this trend. If confirmed, this trend 
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in the rh-L/L Edd relation may naturally explain the above observed lack of objects with very steep 
soft X-ray spectra (i.e. T s > 3.5) at the lowest FWHM end. 



6.3. Comparison with IMBH AGNs 



Three objects in our sample were studied by iMiniutti et al.1 (I2009|) as AGNs with intermedi- 
ate mass black holes (EVIBHs), a term so metimes used to refer to black holes with M B h< 1O 6 M 
in the literature (e.g. iGreene & Hoi 120041) . Since M BH oc FWHM 2 , most IMBH AGNs must have 
broad line widths falling within the conventional criterion of NLSls (< 2000 km s" 1 ), but may not 
necessarily possess the characteristics of typical NLSls, i.e., strong F ell emission, high Edding- 
ton ratios, and significant soft X-ray excess (e.g. Greene & Hoi 2004). The X-ray properties of 



IMBH AGN samples have b e en studied by severa l authors ([Greene & Hoi 12007 al : iDewangan et al. 



20081 : iDesroches et al J l2Q09t, IMiniutti et al.l 120091) . and a large diversity has been found. For in- 
stance, the soft X-ray (0.5-2 keV) photon indices are found to fall into a large range (r s = 1-2.7, 
Desroches et al.ll2009l) . The flat X-ray spectral slopes, as well as some other properties, are very 
similar to those of typical Seyfert galaxies w ith M RH = 10^~ 8 M^ . Some, especially tho se with 
low L/L Edd , do not show soft X-ray excess (|lwasawa et al.1 l2000l : IDewangan et al.l 120081) . as in 
NGC 4395, the prototype of this kind. 

We suggest that IMBH AGNs, albeit their small linewidths as for NLSls, have diverse ob- 
served properties, depending on the Eddington ratio. Those accreting at high L/L Edd values are 
probably more NLSl-like, e.g. the presence of a significant soft X-ray excess, strong Fell, steep 
T s , such as the three IMBHs in Miniutti et al. (2009) and also included in our sample as VNLS1. 
On the other hand, there exists a population accreting at low L/L Edd , which exhibit properties re- 
sembling closely those of classical Seyfert galaxies with more massive M BB , in both optical and 
X-ray (week Fell, relatively flat T s , non-ubiquity of the soft X-ray excess). The observed spectral 
properties of Seyfert galaxies depend much more strongly on mass accretion rate than on black 
hole mass. In this regards, the conventional definition of NLSls may have to be revised. On the 
other hand, NLSls and IMBH AGNs show similar timing property. This is not surprising given the 
postulation that the X-ray variability of AGNs is believed to be largely determined by black hole 
mass, rather than accretion rate. 



7. Summary 



NLSls with very small broad-line widths represent the extreme of Seyfert 1 AGNs, which 
have the largest L/L Edd /M B H ratio among all AGNs known so far. Here we investigated the X-ray 
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properties of a homogeneously selected sample of NLSls with FWHM(H/3) < 1200kms _1 , using 
the archival XMM-Newton data. We note that our sample is not complete in the sense that only 
those observed with XMM-Newton with good spectral S/N ratios are included, which might be 
biased toward relatively bright objects in X-rays. This should be kept in mind when comparing our 
results with the others. 

No significant Fe Ka emission line is detected, which should be at most weak in such objects. 
It is found that the soft X-ray excess is ubiquitous in the objects which have the 0.2-10 keV spectra 
available. The temperatures of this component, when fitted with a blackbody (or disk blackbody) 
model, all fall within 0. 1-0.2 keV, significantly higher than the prediction of the standard disk 
model. Our result highlights the puzzling independence of the thermal temperature on M B h by 
extending it to NLSls with narrower FWHM(H/3), i.e., smaller M BH and/or higher L/L Edd . The 
failure to ascribe the soft X-ray excess to the Wien tail of the disk blackbody emission in these 
VNLSls (with similar M B h and L/L Edd val ues to RX J1633+ 4718, though) highlights the question 
as to why RX J1633+4718 is so unique (see Yuan et aPboiO , for a brief discussion). A range of vi- 



able models, including Comptonization, disk reflection, smeared absorption, and the p-free model 
were used to fit the soft X-ray excess. In general, the disk reflection and Comptonization models 
tend to give the best fits. The relative strength of the soft X-ray excess appears to be independent 
of FWHM(H/3) over a large range, indicating that the excess component is not particularly strong 
in these VNLSls compared to PG quasars with much broader linewidth. 

The soft X-ray spectra in 0.2-2.4 keV have a mean photon index of r s =2.83^o 'o with a large 
intrinsic scatter (cr = 0.3 I^q), while the 2-10keV spectra have a mean T h of 2. 19^ }| with a 
small intrinsic scatter consistent with zero. Thus VNLS 1 s have the spectral slopes in both bands no 
steeper than "normal" NLS 1 s with broader linewidth. There is little or no intrinsic X-ray absorption 
in most of these VNLSls, indicating that the flattening of the r s -FWHM anti-correlation below 



FWHM~1200kms , as suggested in iZhou et al.l (|2006, using a much larger sample but with F 



estimated from hardness ratios), is not caused by absorption but most likely intrinsic. Although 
th is trend is not statis tically significant when combining our current sample with the Swift sample 



of lGrupe et al.l(|2010l) . both with T s derived from spectral fitting, there appears a lack of AGNs with 
both very narrow FWHM(H/3) (< 1000 km s" 1 ) and very steep soft X-ray spectra (i.e., T s > 3.5). 
Furthermore, there is a tentative hint that the hard X-ray slopes r h of our objects fal l systematically 
below the extrapolation of the suggested r h -L/L Edd correlation ( Risaliti et all2009 ) to high L/L Edd 



values. We argue that these two trends, if confirmed, might in fact be driven by the same underlying 
physical process. Similar to other "normal" NLSls, these VNLSls also follow the same a ox -L opt 
relation as for normal Seyferts and quasars. 



All of the sample objects show rapid variability in X-rays, with two-fold timescales of 1-2 
hours. The short variability timescales and the conformance with the variance excess-black hole 
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mass relation for normal Seyferts tends to suggest that the black hole masses in these VNLSls are 
likely truly small, as commonly thought, and the present M BH estimators based on the linewidth- 
luminosity scaling relation is applicable to NLSls. 
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Table 1 . Basic parameters of the sample objects 



No. 


Name 


SDSS Name 


z 


logAL A51 oo 


FWHM(H/3) 


FWHM(Ha) 


F(Ha Ac ) 


R4570 


logMBH 


logL bol /L Edd 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1 


J0107+1408 


J0107 12.0+140845 


0.076 


42.96 


787±31 


709±12 


1500±12 


0.35±0.05 


5.87 


-0.13 


2 


J0740+3118 


J074020.2+311841 


0.295 


44.20 


1135±28 


1090±17 


2121±24 


0.92±0.05 


6.90 


0.073 


3 


J0922+5120 


J092247.0+512038 


0.159 


44.01 


1132±27 


1002±13 


3058±23 


1.37±0.04 


6.63 


0.145 


4 


J0940+0324 


J094057.2+032401 


0.060 


43.12 


1119±95 


810±34 


2249±36 


0.92±0.09 


5.98 


-0.08 


5 


J1000+5536 


J100032.2+553631 


0.215 


43.79 


1065±75 


1216±35 


1056±17 


0.28±0.10 


6.76 


-0.20 


6 


Jl 114+5258 


Jl 11443.7+525834 


0.079 


43.25 


970±60 


866±21 


1158±14 


0.65±0.08 


6.04 


-0.02 


7 


Jl 140+0307 


Jl 14008.7+030711 


0.081 


43.12 


675±41 


571±18 


1337±17 


0.97±0.08 


5.69 


0.200 


8 


J1231+1051 


J123126.5+105111 


0.304 


43.92 


957±23 


1200±332 


1143±255 


0.57±0.05 


6.93 


-0.23 


9 


J 1246+0222 


J124635.2+022209 


0.048 


43.49 


811±27 


709±15 


11669±118 


0.83±0.04 


6.09 


0.171 


10 


J1331-0152 


J133 141. 0-015213 


0.145 


43.56 


1192±42 


1044±16 


1929±15 


0.35±0.03 


6.54 


-0.20 


1 1 


J1357+6525 


J135724.5+652506 


0.106 


43.14 


737±41 


694±16 


1471±16 


0.45±0.07 


6.00 


-0.08 


12 


J1415-0030 


J141519.5-003022 


0.134 


43.36 


1045±27 


954±11 


1735±13 


0.76±0.04 


6.40 


-0.26 


13 


J2219+1207 


J221918.5+120753 


0.081 


43.66 


982±38 


886±15 


4369±38 


1.11 ±0.06 


6.32 


0.115 



Note. — Col. (2): abbreviated name of objects; Col. (3): SDSS name; Col. (4): redshift; Col. (5): monochromatic luminosity at 5100 A (ergss -1 ); Col. (6): H/3 
linewidth (km s~'); Col. (7): Ho linewidth (km s~'); Col. (8): Ha broad component flux (10~ 17 ergs s~' cm -2 ); Col. (9): the optical Fe II strength relative to the H/3 
broad component; Col. (10): black hole mass (Mg); Col. (11): Eddington ratio 
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Table 2. Log of XMM-Newton observations 



Name 


Date 


Off-Axis 




Exposure time 




Note 








PN 


MOS1 


MOS2 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


J0107+1408 


2005-07-22 


0.0 


15.1 


27.1 


27.1 


c 


T0740+3 1 1 8 


9001 -04- 1Q 


1? 7 

1 _■ / 








o 


J0922+5120 


2005-10-08 


0.0 


5.2 


19.2 




c 


J0940+0324 


2005-10-30 


11.8 


21.9 


26.1 


26.2 


b 


J1000+5536 


2001-04-13 


12.8 






8.2 


a 




2003-10-14 


13.5 


14.2 






a 


Jl 114+5258 


2003-04-25 


11.3 


3.9 


6.6 


7.0 


a 


Jl 140+0307 


2005-12-03 


0.0 


30.8 


38.5 


39.2 


c 


J1231+1051 


2003-07-13 


12.0 




45.3 




a 




2005-12-13 


12.0 






68.2 


a 




2005-12-17 


12.0 




91.6 


91.6 


a 


J1246+0222 


2001-06-17 


0.0 


3.1 






c 


J1331-0152 


2001-07-29 


12.0 




32.3 


32.3 


b 


J1357+6525 


2005-04-04 


0.0 


14.5 


21.1 


20.6 


c 


J1415-0030 


2003-02-08 


10.0 


9.5 


13.7 


14.2 


b 


J2219+1207 


2001-06-07 


0.1 


7.2 






c 



Note. — Col. (2): observation date; Col.(3): off axis angle in arcmin; Col.(4)-(6): 
cleaned exposure time of the three EPIC cameras in kilo-second; Col. (7): energy 
range in which source spectrum is extracted, (a) 0.2-2.4 keV, (b) 0.2-7 keV, (c) 
0.2-10 keV. 
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Table 3: XMM-Newton spectral fits with an absorbed power-law model 



Name 
(1) 


(2) 


N£ 
(3) 


r s0 ft 

(4) 


log Ljoft 
(5) 


xi/dof 
(6) 


rWa 

(V) 


log Lhard 

(8) 


(9) 


EW 
(10) 


J0107+1408 


3.37 


4.2 ±0.8 


2.63 ±0.08 


43.07 


1.0/230 


2.26±0.22 


42.47 


1.3/44 


<45 


J0740+3118 


4.32 




3.38 ±0.30 


44.72 


c-stat 










J0922+5120 


1.20 


2.7±0.5 


3.72±0.07 


44.23 


1.7/226 


2.28±0.36 


43.31 


1.3/24 


<228 


J0940+0324 


3.25 




2.43 ±0.03 


43.00 


1.0/244 


20 2+o.i9 


42.56 


1.1/39 


<497 


J1000+5536 


0.83 




2.03±0.38 


43.33 


c-stat 














2.34 ±0.25 


43.14 


1.5/18 










Jl 114+5258 


0.99 




2.77 ±0.07 


43.16 


0.8/60 










Jl 140+0307 


1.88 


1.3±0.2 


2.89 ±0.03 


43.44 


1.3/375 


2.06±0.13 


42.65 


1.0/94 


< 1 1 1 


J1231+1051 


2.40 




2.89±0.18 


44.10 


0.9/26 
















2.80±0.16 


43.83 


1.5/28 
















2.90 ±0.07 


43.93 


1.0/69 










J 1246+0222 


1.74 


1.83±0.5 


2.97 ±0.06 


43.65 


1.2/198 


9 oc+0.35 
z "^ J -0.23 


42.74 


1.1/17 


<214 


J1331-0152 


2.33 




2.67 ±0.05 


43.65 


1.0/58 


1 nc+0.69 
ly;, -0.35 


43.18 


0.8/11 


<361 


J 1357+6525 


1.21 


2.3 ±0.4 


2.68 ±0.06 


43.43 


1.1/265 


O 90+O.I6 

y -0.32 


42.85 


1.1/46 


<276 


J1415-0030 


2.84 




2.76 ±0.06 


43.49 


1.0/104 


9 91+0.51 
z - z '-0.47 


42.81 


1.4/11 


<842 


J2219+1207 


4.81 


2.8±0.7 


3.15 ±0.07 


44.19 


1.0/264 


2.39±0.16 


43.21 


1.1/38 


<125 



Note. — Col. (2): Galactic column density in 10 20 cm" 2 ; Col. (3): column density of intrinsic neutral absorption 
in the object's rest frame in 10 20 cm" 2 ; Col. (4): fitted power-law photon index in the soft X-ray band (0.2-2.4 keV); 
Col. (5): absorption corrected luminosity in 0.2-2.4keV in ergss" 1 ; Col. (6): reduced x 2 ; Col. (7): fitted power-law 
photon index in the hard X-ray band (mostly in 2-10 keV; see Tabled; Col. (8): absorption corrected luminosity in 
2-10keV in ergss" 1 ; Col. (9): reduced x 2 ; Col.(10): rest frame equivalent width of the Fe Ka line in units of eV 
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Table 4:: Results of spectral fits in 0.2-10keV with different 
models for the soft excess 



Blackbody 



Name 


r 


kT(keV) 


BB/total (a) 


X 2 /dof 






JO 107+ 1408 


9 og+0.08 


0.11 ±0.02 


0.166 


277/264 






J0922+5120 


3.18±0.09 


0.12±0.01 


0.380 


356/255 






J0940+0324 


2.29±0.05 


1 1 +004 

u - 1 i -0.03 


0.079 


278/273 






Jl 140+0307 


2.56±0.04 


0.14±0.01 


0.206 


570/472 






J1246+0222 


2.61±0.02 


0.15±0.01 


0.253 


220/212 






J1331-0152 


2.27±0.06 


0.09±0.01 


0.186 


61/68 






J1357+6525 


2.36±0.09 


0.14±0.01 


0.219 


319/303 






J 14 15-0030 




0.10±0.02 


0.160 


102/109 






J2219+1207 


2.70±0.04 


0.14±0.01 


0.208 


293/293 






Comptonization 


Name 


r 


W plasma (keV) 


T 


X 2 /dof 






JO 107+ 1408 


2.32±0.06 


0.25±0.07 


10.9±4.3 


272/263 






J0922+5120 


2.09±0.19 


0.17±0.01 


1 Q Q+25 


325/254 






J0940+0324 


2.04±0.19 


0.30±0.10 


13.5±2.1 


269/272 






Jl 140+0307 


2.19±0.15 


0.22±0.03 


18.7±7.9 


540/471 






J 1246+0222 


2.34±0.23 


0.21±0.02 


21.0±5.5 


217/211 






J1331-0152 


2.12±0.21 


0.17±0.07 


17.1±7.9 


60/67 






J1357+6525 


2.22±0.04 


0.19±0.02 


23.5±7.0 


320/302 






J 14 15-0030 


2.04±0.20 


0.20±0.06 


20.0±4.8 


99/108 






J2219+1207 


2.36±0.15 


0.21±0.04 


18.1±7.0 


279/292 






Disk Reflection 


Name 


r 




log£ 


kT (keV) (fo) 


Flux Frac. (c) 


X 2 /dof 


JO 107+ 1408 


2.24±0.17 


4.93+3.31 


3.45±0.4 


0.05±0.008 


0.3 


260/259 


J0922+5120 


2.41±0.25 


U3+11.07 


o OO+0.13 
J - JJ -0.12 


0.06±0.01 


0.5 


257/250 


J0940+0324 


2.16±0.12 


1 41+8.98 


3J6\ 2l 




0.9 


273/271 


J 1140+0307 


2.38±0.08 


3.86! 089 






0.9 


524/470 


J 1246+0222 


2.43±0.10 


3.82! 232 


3.75_o 19 




1.0 


212/210 


J1331-0152 


2.29±0.25 


1.24+ 315 


3-27^25 




0.3 


69/66 


J1357+6525 


2.19±0.07 


6 05 +L91 


o to+0.30 
■ 3 - zo -0.23 




0.5 


320/301 


J 14 15-0030 


2.24±0.24 


1.25! 2 ' 24 


3 i O+0.12 




0.7 


104/109 


J2219+1207 


2.43±0.05 


4.46+2.15 


3-79\ 15 


0.08^ 


0.8 


272/289 



Smeared Absorption 
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Table 4:: 


(continued) 






Name 


r 


log£ N H 


(10 22 cm" 2 ) 


a (v/c) 


X 2 /dof 


JO 107+ 1408 


2.47±0.21 


3.00±0.71 




0-50V 6 


278/264 


J0922+5120 


2.69±0.07 


3.11 ±0.04 


50! 2 


5 + 

u - J -0.05 


388/255 


J0940+0324 


2.22±0.10 


3.46±0.35 


35^4 


0-44! . 28 


275/273 


Jl 140+0307 


2.37±0.03 


3.40±0.05 


50! 3 


0-5!*!o.05 


607/472 


J 1246+0222 


2.56±0.09 


3.58±0.10 


50! 8 


47 + 


291/211 


J1331-0152 


2.42±0.07 


2.75±0.36 


14! 25 


5+ 


60/67 


J1357+6525 


2.26±0.06 


3.59±0.06 


50 + _ 9 


5+ 


359/303 


J 14 15-0030 


2.60±0.12 


3.52±0.50 


9 + - 5 


0-49!o.i6 


112/110 


J2219+1207 


2.60±0.20 


3.43±0.03 


50U 


5+ 

u - J -0.05 


320/292 



Notes. For SDSS J0107+1408 with free absorption in the fitting. For SDSS J1415-0030 

an additional ionized absorption is applied (see text). The blank parameter errors denote that the 
upper or lower limits are outside of the tabulated parameters range, which are considered to be not 
physically meaningful. (a) Luminosity ratio of the blackbody component to the total component 
in 0.5-2 keV. (fc) Inferred temperatures of an additional blackbody component. (c) Flux ratio of the 
reflected component to the total component in 0.2-10 keV. 
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0123 56789 -2 -1012 

R4570 'og M BH (M ) log (L bol /L E<J d) 



Fig. 1. — Distributions of the Fe II to H/3 flux ratio (R4570), black hole mass and the Eddington 
ratio for our sample (red), the total Zhou'06 NLSls sample (black) and the whole subsample with 
FWHM(H/?) < 1200 km s" 1 (blue). 




Fig. 2. — Composite optical spectrum of our sample objects derived from their individual SDSS 
spectra. The strong Fe II multiplet emission is characteristic of typical NLS 1 spectra. 
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0.8 [ 




<r> 

Fig. 3. — Confidence contours (at the 68%, 90% and 99% confidence levels) of the mean and stan- 
dard deviation of the intrinsic distributions (assumed to be Gaussian) of the soft (r s ) and hard (1^) 
X-ray photon indices for our VNLS1, which are derived using the Maximum-likelihood method 
(see text). Pluses indicate the best-estimated values. 
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Fig. 4. — XMM-Newton spectra of the four objects among our VNLS1 sample (PN spectra except 
for J133 1 -0152 of which MOS spectra are used). The power-law model fitted to 2-10keV spectra 
and its extrapolation to the soft X-ray band is also shown. Soft X-ray excess emission is clearly 
present. 
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Fig. 5. — The observed temperature of the soft excess is plot ted versus the black ho le mass. Filled 
circles are our res ults, stars for the rad io-quiet PG quasars (|Piconcelli et al.ll2005|) and diamonds 
for type 1 AGNs ([Crummy et al.1 120061) . The dotted-dashed lines are the maximum temperature 
expected from the accretion disc. 
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Fig. 6. — Soft X-ray excess strength, parameterized as the ratio of the blackbody to the total lumi- 
nosity in the 0.5-2 keV range, versus the lin ewidth for the VNLSls in our sample (open circles) 
and the radio-quiet PG quasars (diamonds) in lPiconcelli et al.l (|2005l) . 
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Fig. 7. — Residuals of spectral fits to the XMM-Newton spectra with various models to account for 
the soft X-ray excess, which are, from top to bottom for each panel, blackbody, Comptonization, 
disk reflection, and smeared absorption model. An additional blackbody component is added in 
the disk reflection model for JO 107+ 1408, J0922+5120 and J2219+1207. 
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Fig. 8. — Soft X-ray pho ton index r s versus H/3 linewidth relation for our VNLSls (filled circles). 
Stars are the results from lGrupe et al.l (|2010t ). 
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Fig. 9. — XMM-Newton X-ray lightcurves in the 0.2-10keV band for five of the VNLSls in our 
sample. The time binsize is 200 s, except for SDSS J133141.0-015213 (500 s). 
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Fig. 10. — Dependence of a ox on the 2500 A monochromatic luminosity for the VNLSls in our 
sample (filled dots) and the ordinary NLSls (open squares) in the sa mple of iGallol (I2006ah . The 
solid line represents the relation for radio-quiet type 1 AGNs given in lStrateva et al.1 (|2005l) . 
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Fig. 11. — Relationship between the hard X-ray (2-10 keV) photon index and the Eddington ratio 
for our sample objects. Filled circles represent the results from the fits with a simple power- 
law model and open circles from a pow er-law plus disk refl ection model. The solid line is the 
extrapolation of the relation suggested by iRisaliti et al.1 (I2009f) (M B h estimated from H/3, the same 
as in our paper) and the dash-dotted line represents the dispersion. The vertical dotted line marks 
the higher end of the L/L Edd range in the Risaliti's sample. 
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A. COMPARISONS WITH PREVIOUS RESULTS FOR INDIVIDUAL OBJECTS 



For six objects in our sample, the XMM-Newton data have been presented previously. Here 
we compare our results with those in the literatures. 

J0107+1408 , J140 +03 07, J1357+6525 — The se three objects were studied as IMBH AGNs 



by iMiniutti et al.1 (120091) and iDewangan et al.l (120081) . For the model fits with blackbody, ionized 
disk reflection, and smeared absorption, our results are consistent with the previous ones e xcept that 
the disk ionization parameters of ours are somewhat higher than those in IMiniutti et al.1 (|2009l) . In 
addition, we fit the spectra with the Comptonization and 'p-free' model, which was not considered 
in the previous papers. 

J 1246+0222 — The fitting re sults for the soft X- r ay excess componen t are consistent with 
those presented in the literatures (IPorquet et al.l 120041 : ICrummy et al.l 120061 : iMiddleton & Done 
20071). However, in the sim ple power-law fitting we find a soft X-ray photon index of 2.97±0.06 
while IPorquet et all hoo4) gave the value of 3.72±0.09. 



J1415-0030 — The XMM-Newton data was presented previously by iFoschini et al.l (|2004l) . 
who fitted the spectrum with the simple power-law plus blackbody model only, and gave a some- 
what flatter photon index (1 .8±0.2) than ours (2.42^ jq). We find that there are an edge-like feature 
around 0.6 keV for J141 5-0030. 



J2219+1207 — The spectrum was fitted with a power-law plus blackbody model by lGallo et al 



d2006bl) . giving a result consistent with ours. The authors also noted the broad excess emission 
feature around 5.8 keV, which was fitted with a broad emission line; while in our work the disk 
reflection model can well reproduce this feature. 



